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Abstract 

Recently strong electric fields (up to 10 9 V/cm) have been discovered, which affect 
the neutrons moving in noncentrosymmetric crystals. Such fields allow new polarization 
phenomena in neutron diffraction and optics and provide, for instance, a new feasibility 
of a search for the neutron electric dipole moment (EDM). A strong interplanar electric 
field of the crystal and a sufficiently long time for the neutron passage through the 
crystal for Bragg angle close to tt/2 in the case of Laue diffraction make it possible to 
reach the sensitivity achieved with the most sensitive now magnetic resonance method 
using ultra cold neutrons (UCN method). 

A series of experiments was carried out in a few last years on study of the dynamical 
diffraction of polarized neutrons in thick (1-10 cm) quartz crystals, using the forward 
diffraction beam and Bragg angles close to 90° . As well new neutron optics phenomena 
were investigated. The following effects was first observed: 

• the effect of an essential time delay of a diffracting neutron inside the crystal for 
Bragg angles close to it/2; 

• the phenomenon of a neutron beam depolarization in time of Laue diffraction by 
noncentrosymmetric Q-quartz crystal; 

• the effect of a neutron-optical spin rotation for neutrons moving through a non- 
centrosymmetric crystal with the energies and directions far from the Bragg ones; 

• the controlled spin rotation effect for the Bragg reflected (from a smal part of 
crystal near its exit face) neutrons in a slightly deformed quartz crystal. 

The feasibility of experiment on a search for neutron EDM using Laue diffraction in 
crystals without a center of symmetry was tested at the reactors: WWR-M in Gatchina 
and HFR in Grenoble. It was shown that the sensitivity can reach (3 — 6) • 10 -25 e • cm 
per day for the available quartz crystal and cold neutron beam flux. 



1 Introduction 

The origin of CP-symmetry violation (where C is a charge conjugation and P is a spatial 
inversion) is of a great interest since its discovery in the decay of neutral K-mesons about 
40 years ago. CP violation leads in turn to the violation of the time reversal symmetry (T) 
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through the CPT invariance (CPT-theorem). Existence of nonzero neutron EDM requires 
violation of both P and T invariance. Different theories of CP violation give widely varying 
predictions for a neutron EDM (see for instance 012)- Therefore, the new experimental 
limits on the EDM value would be of great importance for understanding the nature of the 
CP violation as well as of the Universe baryon asymmetry, which are beyond the Standard 
Model. 

The most precision method of the EDM measurement now is the magnetic resonance 
one PIE El 1 using the ultracold neutrons (UCN method). 

Here we shall discuss another possible way for the neutron EDM search, using Laue 
diffraction of the cold neutrons in the noncentrosymmctric crystal [B]. Earlier we have 
shown that the neutron, moving through a noncentrosymmetric crystal, may be influenced 
with a strong intcrplanar electric field. Its value depends on direction and value of a neutron 
wave vector, reaching a maximum (up to 10 9 V/cm), when the Bragg condition for some sys- 
tem of crystallographic planes is satisfied. These fields result in new polarization phenomena 
observable in neutron diffraction and optics. In particular, the Schwinger interaction of the 
neutron magnetic moment with such fields leads to spin dependence of the neutron Pen- 
dellosung picture [7] . The Pendellosung phase shift accompanying the spin flip was measured 
first in ref. [7], and so the interplanar electric field was determined in the experiment on 
dynamical Laue diffraction of polarized neutrons for the (IlO)-plane of the a-quartz crystal. 
Experimental value of the field was obtained to be equal i?mo) = (2.10 ±0.12) x 10 8 V/cm 
and had coincided with the theoretical one. 

The interplanar electric fields are more than four orders of magnitude higher than those 
used in the UCN methods H El of the neutron EDM search. So it was a natural idea 
(BUSED (arisen anew after Shull and Nathans ^5) to use these crystal fields for searching 
the neutron EDM (see also |12l However, the value of the crystal field turned out to 

be still insufficient to reach the sensitivity of the UCN method, but it was shown [Hj that 
use of Laue diffraction for EDM measurements with the Bragg angles close to 90° would 
essentially increase the time, the neutron spends in crystal under the strong electric field. 

Two variants of the method, using the Laue diffraction of polarized neutrons in the 
crystals without a centre of symmetry, were proposed for this purpose. One of them is the 
double crystal variant jS] based on a spin dependence of the Pendellosung picture phase 
and the second is the polarization method [5] using the depolarization effect for neutrons 
diffracting in the noncentrosymmetric crystals. It was shown EI that the sensitivity 
of the method to measure the neutron EDM may be increased more than by an order of 
magnitude by a choice of Bragg angles close to 7r/2, so it may exceed (with the higher 
luminosity taken into account) the sensitivity of the UCN method EI El • That is possible 
for Laue diffraction scheme only. But only the experimental study of these effects can answer 
the question on actual sensitivity of the method to measure the neutron EDM. 

A hypothetical idea to use the interplanar electric fields (if they did exist) for the neutron 
EDM search was discussed in the review |13j . But such crystal properties were not known 
at that time. 

The importance to consider the crystal noncentrosymmetricity was previously pointed 
out in ref. Authors were the first who paid attention to an existence of the 

J Last experimental result for a neutron EDM obtained by this method < 6.3 ■ 10~ 26 e-cm at the 

90% confidence level 
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interference between the electric and nuclear structure amplitudes for neutrons diffracted 
by noncentrosymmetric crystal and proposed to study the Schwinger interaction using such 
crystals. 

In the work |1U) the effect of neutron spin rotation due to such interference has been 
discussed. Similar, but more detailed theory of neutron optical activity and dichroism for 
diffraction in noncentrosymmetric crystals has been developed in ref. |f 2) . 

It has been shown ^U] that the spin rotation effect in a non-absorbing crystal can take 
place only for Bragg scheme of diffraction. The deviation from the Bragg condition by about 
the Bragg width is necessary in this case to observe the effect for the transmitted beam. 
That reduces the effect, because it is proportional to l/vl + w 2 , where w is the parameter 
of angular deviation (measured in the units of Bragg halfwidth) from the Bragg direction, 
see 8 . The nuclear absorption is necessary to have a spin rotation effect in the case of Laue 
diffraction QUIU2]. 

The possibility of a search for neutron EDM by measuring a spin rotation angle was 
analyzed for the Bragg diffraction scheme ^0] and for the Laue scheme |12| . 

We have shown jjTf that for some system of crystallographic planes in noncentrosymmetric 
crystals the positions of the electric potential maxima can be shifted relative to those of a 
nuclear potential. Hence, the diffracted neutrons will move in the crystal under a strong 
(10 8 — 10 9 V/cm) interplanar electric field (because of the neutron concentration on the 
nuclear potential maxima or between them) . Such a concept turned out to be very fruitful for 
further consideration and understanding the different phenomena, concerning the neutron 
diffraction and the optics in noncentrosymmetric crystals |S| . For example it allowed to 
predict and to give a simple description of such new effects as the spin dependence of the 
pendulum phase, the depolarization of the diffracting neutron beams, the independence of 
the effects due to the Schwinger interaction on the neutron wavelength and Bragg angle for 
given crystallographic planes etc. p iT^llTH] . 

We should note also, that the Bragg reflection of neutrons from the neutron-absorbing 
centrosymmetric crystal of CdS was used earlier for EDM search. But the sensitivity 
of this method is much lower than that of UCN-method, because the depth of neutron 
penetration into the crystal, which determines a time the neutron stay in crystal, was very 
small (about 7T0~ 2 cm). Now a new variant of the method is proposed and developed, 
using the multiple reflections from the silicon centrosymmetric crystal |f 71 II8| . 

In the work it was reported that the effect of the neutron spin rotation has been ob- 
served due to the spin-orbit (Schwinger) interaction, using the Bragg scheme of the diffrac- 
tion in the noncentrosymmetric crystal with a small deviation of a neutron momentum 
direction (by about a few Bragg width) from the Bragg one, because the effect disappears 
for the exact Bragg direction in this case. However, the experimental value of the spin rota- 
tion angle |19j turned out to be a few times less than the theoretical one. Authors were in 
difficulty to explain the origin of such discrepancy, but it was very likely due to imperfection 
of the used crystal. 

Here we consider neutron-optic effects for neutrons, moving through a noncentrosymmet- 
ric crystal with the energies and the directions far from the Bragg ones, when the deviation 
from the exact Bragg condition reaches (10 3 — 10 5 ) Bragg widths. 

The theoretical estimations have shown |20j. that for the polar noncentrosymmetric 
crystal (PbTiO^, for instance) the value of a resultant electric field, acting on a neutron, 
can reach ps 2 x I0 6 V/cm for wide range (about four orders exceeding the Bragg width) 
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of the neutron directions and wavelengths. Such a field is a result of superposition of the 
fields from a few different crystallographic planes. Therefore, the spin rotation effects may 
be not too small in neutron optics in comparison with the diffraction case. 

The observation of such effects may be of interest for development of new methods for 
searching the neutron EDM. It also may be useful for experimental searching the T-odd 
part of the nuclear interaction, using neutrons with energy close to the P-resonance one 
|21|. because it is hardly possible to observe dynamical diffraction in a large crystal for 
neutrons with the energy ~ 1 eV because of too high requirements to a crystal quality. 

The phenomenon can be used also for measurements of the interplanar electric fields 
affecting the neutrons in the crystals without centre of symmetry A new specific spin 
neutronography arises in this case for crystallographic planes with nonzero electric fields. 



2 Diffraction in a noncentrosymmetric crystal 

As it follows from the dynamical diffraction theory, a movement of a neutron through the 
crystal in a direction close to the Bragg one for some system of crystallographic planes can be 
described by two kinds of Bloch waves ip^ and (see for example, formed as a result 
of neutron interaction with the periodic nuclear potential |S] V^r) = Vjj^ + 2v^ cos(gr), 
where g is a reciprocal lattice vector describing the system of crystallographic planes, |g| = 
2ir/d, d is an interplanar spacing, Vq is the average nuclear potential of the crystal, 

= cos7e 4fc<1>r + sin 7 e lfe l 1>r , (1) 



V> (2) =-sm 7 e tk r +co S je lk y r , (2) 

where tg2 7 = \U g \/A = 1/w = 7B /ft; < 7 < tt/2, = 2mv g v /h 2 , fe^' 2) = fc (1 ' 2) + g; 
A = (kg — k 2 )/2 . The parameter w is the ratio of angular deviation from the exact Bragg 
direction fl = 9 — 9 b to the angular Bragg width 75 , therefore w describes the relative 
deviation from the exact Bragg condition. 

The intensities of direct and reflected waves in the states QJ, J5J are given by 

9 l r A a 1. 

cos 2 7 =-[l+ - 9 ] = -[1 



The expressions JQ, describe two standing waves (in the direction g normal to the 

crystallographic planes), which are moving in the direction kfi 1 = k^ 1 ' 2 ) + g/2 along the 

planes (see Fig. QJ. A small difference of the wave vectors k^ x \ k^ is a result of neutron 
concentration on "nuclear" planes (for and between them (for ip^), 



kW 2 = K 2 -A±JA 2 + \U g \ 2 . (3) 



where K 2 = k 2 + = 2m(E + V^)/h 2 , m, E, k are respectively the mass, energy and 
wave vector of the incident neutron. 
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Figure 1: Passage neutron through the crystal. 

The values Vq* , Vg have an order of 10 -7 eV, so for thermal and cold neutrons with 
the energies of (10 -1 — 10~ 3 ) eV we can consider that k^ w w k. The propagation 
velocities of these waves along the crystallographic planes are 2 

J 1 * 2 ) = l\ k M +g/2 | = l k (^) cos e B ^vcos9 B , (4) 
11 m m 

where v = fikjm — 2-KTi/(\m) — tt h/{md sin 9 b) is the velocity of the incident neutron, 9b 
is the Bragg angle, A = 27r/fc is the wavelength of the incident neutron (A = 2d sin 9b)- A 
number of the dynamical diffraction phenomena (see, for example, |22l 1231 123] '). including 
effects caused by the neutron EDM (HUH], are determined not by a total neutron velocity v, 
but its component along the crystallographic planes Vu = v cos In particular, the time the 
diffracting neutron spends in crystal, which equal to tl = L/(vcos9b) ~ L/[v(n/2 — 9b)}, 
sharply grows for Bragg angles close to 7r/2, where L is the thickness of a crystal. That 
allows to increase the sensitivity of the diffraction method to neutron EDM at least by an 
order |U] . Therefore the values 3 Et can be of the same order for UCN and diffraction 
method (for Bragg angles sufficiently close to 7r/2) [HUHj despite the fact that the storage 
time for UCN (~ 100 s 0) is essentially more than the time of a neutron passage through 
the crystal 4 . 

We proposed the polarization method for searching a neutron EDM based on the pre- 
dicted effect of depolarization of the diffracted neutron beam [5] in noncentrosymmetric 
crystal. This method has some advantages (such as its relative simplicity and less sensitiv- 
ity to crystal imperfection) over the method based on the spin dependence of the neutron 

2 Here we neglect the Pendellosung oscillations (arising from the interference of waves of different type) 
because in our case they are averaged over Bragg angles for a slightly divergent beam. 

3 The set-up sensitivity is determined by u(D) oc 1/Ery/N, where cr(D) is an absolute error of EDM 
measurement, N is a total value of accumulated events. 

4 The ways of possible improving the UCN-method are discussed |25l . 
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Pendellosung phase [Sj. 

Essence of the effect is as follows. In the noncentrosymmetric crystal the diffracting 
neutrons in two Bloch states are moving under opposite electric fields [7| 0D because of 
the shift of the "electric" planes relative to the "nuclear" ones, see Fig. ^ We mean that 
" nuclear" or " electric" planes are determined by the maximum positions of the corresponding 
periodic potentials of this plane system 

V g N (r)=V Q N +2v? cos(gr), 

V g E (r)=V E + 2vf cos(gr + 0f), 

the shift <j) E may be calculated using the crystal structure. The electric field of the plane 
system has a form 

E(r) = -grad V g E (r) = 2vf gsin(gr + <f>f). 
So the mean electric fields acting on a neutron in the states ijj^ and ip^ will be equal 

<V,«|E(r)|^>) = -(^ 2 )|E(r)|V^} = -j^L= = g f . (5) 

where 

E s = g*fsin0f, (6) 

is the maximum electric field acting on a neutron under the exact Bragg condition (w = 0). 

Therefore spins in these states will rotate in the opposite directions due to Schwinger 
interaction, that in turn will lead to a decrease of the neutron beam polarization (see Fig. |2J ■ 
If an initial spin orientation is normal to the "Schwinger" magnetic field = [E fl x V||]/c, 
then for the case w « 1 (exact Bragg condition) the spin rotation angle in both states will 
be equal to PITS] 



^H s g L eEgL 

-fc — ^Mn o 

nvn m v c z 



A^o = = ±Mn-^r, (7) 



because E g _L vy and H g = E g v\\/c. Here the signs ± are related to different states (1), (2) 
respectively, fj, n = —1-9 is the neutron magnetic moment in nuclear magnetons. As a result 
the value of neutron beam polarization P will depend on A0q in the following way 5 : 

P = P cos = P cos f^M] , ( 8 ) 

Pq is the incident beam polarization (see Fig.[5J). 

The polarization P can be decreased down to zero by a choice of such a crystal thickness 
Lq that makes the spin rotation angles equal to ±7r/2. Theoretical calculation for (110)- 
planes of a-quartz gives Lq — 3.5 cm. 

The existence of the neutron EDM leads to a slight polarization Pedm along , equal to 

4DE g Lo 4D c I 

Pedm = r- 2 — = 5- oc — — — , (9) 

irnvt\ ir[i vcosob tt/2 — Ub 



5 This result is obtained averaging the Pendellosung oscillations over Bragg angles. The angular period 
of this oscillations in our case is ~ 10 — 5 rad and the angular divergence of the neutron beam is ~ 10 -2 rad. 
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Figure 2: The neutron spin for two Bloch waves (shown schematically as the red and 
blue circles) will rotate in opposite direction under opposite fields. When the angles of 
spin rotation become equal to 7r/2, the both diffracted neutron beams will be depolarized 
entirely. The existence of neutron EDM will result in a polarization along the Schwinger 
magnetic field. The sign of this polarization will be different for two crystal position A and 
B. 
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because cos9b oc tt/2 — 9b for 9b — > 7r/2. Here D is the neutron EDM. The turn of the 
crystal by the angle 29b (see Fig. 01 will change the Pedm sign but will not do that for 
residual polarization. High precision of the crystal turning (better than 10 -5 rad) gives the 
possibility to exclude the systematic errors and to select the EDM effect. 

The principal scheme of the Laue diffraction method to search for the neutron EDM 
|5ni|^7] is shown in Fig. EI For two crystal positions R and L with the same Bragg angle but 
with opposite directions of the electric held, the polarization Pedm will have opposite signs 
whereas a residual polarization will have the same sign for both crystal positions. Therefore 
we should put the initial neutron spin along the neutron velocity (Y axis) and compare the 
components of the polarization along the Z axis for two crystal position marked by R and 
L in Fig. El 



Figure 3: Scheme of the experiment for a neutron EDM search by Laue diffraction. The 
presence of a neutron EDM will lead to a small Y-component of the polarization, which will 
have different signs for the two crystal positions R and L. 

As it follows from (0 the effect due to the Schwinger interaction does not depend on 
such neutron properties as the energy, wavelength and the Bragg angle. It is determined 
by the property of crystal and by the fundamental constants only. For given crystal it is 
the same for any Bragg angles. That gives an additional way to eliminate a false effect 
concerned with the Schwinger interaction by carrying out the measurements for two Bragg 
angles, for example. The EDM effect (0 (in contrast to the Schwinger one) depends on the 
Bragg angle. It essentially increases for 9b — > 7r/2. 

3 Measurement of time the neutron spends in crystal 

Two experiments described below on observing the effect of the neutron time delay in the 
crystal and the depolarization effect were carried out at the WWR-M reactor in PNPI 
(Gatchina, Russia). The scheme of experimental setup is shown in Fig.Qf) [281 I29| . 

The neutron beam formed by neutron guides 1,2 is diffracted by the noncentrosymmetric 
a-quartz crystal 6 (the reflecting (110) planes are normal to the large crystal surfaces) and 
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is registered by the detectors 11. 




Figure 4: Two modifications of the experimental set-up: a) for time-of-fiight measurements, 
b) for measurement of the depolarization effect. 1 is a reactor interchannel neutron guide, 
2 is a multislit polarizing neutron guide, 3 is the BeO polycrystal filter (120 mm), 4,9 
are spin-guide coils, 5,8 are spin-rotation coils, 6 is the a-quartz single crystal, (sizes are 
14 x 14 x 3.5 cm 3 ), 10 is a double multislit polarizing neutron guide, 11 are neutron detectors, 
12 is a beam chopper. A and B are two crystal positions with the same Bragg angles, g is 
the reciprocal lattice vector for the (llO)-plane, Hl is the guiding magnetic field, I is the 
TOF length. 

All neutrons diffracted by the different crystallographic plane systems (for which the 
Bragg conditions are satisfied) give the contribution to the intensity of the direct diffracted 
beam. We used the time-of-flight (TOF) technique to select the specified reflection. The 
mechanical beam chopper 12 was placed before the crystal. The typical TOF spectrum is 
shown in Fig. |SJ The peaks corresponding to neutrons diffracted by the different crystallo- 
graphic planes are clearly visible in the figure. 

If the crystal is located between the beam chopper and the neutron detector, the total 
time of flight of neutron with the wavelength A = 2d sin 9 b will be: 

/ L dm ,, „ „ , . 

r f = T l +T L = - + — = —(lsm8 B + Ltg0 B ), (10) 

v v cos Ob nit 

where t; is the neutron time of flight for a distance I, tl is the time a neutron spends in 
the crystal, L is the thickness of the crystal, 9b is the Bragg angle (for the (110) plane of 



9 




Figure 5: TOF spectrum of the forward diffracted neutrons for Bragg angle 9b = 75°. n T is 
the order number of the TOF channel. The width of the TOF channel is equal ~ 51.2 /is. 
N is the number of accumulated events. 



the a-quartz crystal d — 2.4564A). As follows from (|10|l the time tl of a neutron delay in 
the crystal depends on the Bragg angle as tgOs, while the t; oc sin6*s, so T£ may give a 
very essential contribution to the total TOF of neutrons t/ for 9b close to 7r/2, because 
TL /Ti~L/[l(Tr/2-e B )}. 

A dependence on the Bragg angle of the neutron TOF for forward beam diffracted by 
the (llO)-plane is given in Fig.EJ 

A solid curve describes the calculated dependence Tf{9 B ) (see iJIOJ), and the dotted curve 
does that for ti(0b)- One can see a good agreement between experimental (black points) 
and theoretical dependencies 7y(0s). A control experiment, when the chopper of a neutron 
beam was placed between the crystal and detector, was also carried out. In this case the 
delay of neutron in the crystal does not give a contribution to the measured value and the 
position of the line ti(9b) for (llO)-reflection should coincide with the dotted curve. That is 
what we have observed experimentally (open points). On insertion in Fig. [S]the theoretical 
and experimental dependencies tl(9b) are shown. 

So the experiment have proved that the time of neutron delay in the crystal is not 
determined by the total neutron velocity v, but its component Vn along the crystallographic 
plane. In particular, for 9b = 87° we find that = (0.90 ± 0.02) ms and Um = (39 ± 1) m/s, 
while v — 808 m/s. 
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Figure 6: Dependence of the neutron TOF on the Bragg angle for forward diffracted beam. 



4 Measurement of the depolarization effect 

So we have shown that the polarization of the diffracted beams will be zero for crystal 
thickness L = 3.5 cm, if the initial neutron polarization is perpendicular to the Schwinger 
magnetic field , while for initial polarization parallel to will stay completely polarized. 
We used a set of different angles between directions of initial neutron polarization and 
Schwinger magnetic field to study this effect in detail. 

The scheme of the experimental setup is shown in Fig.0f\ see also [551 EH- The polar- 
ization vector Po of a neutron after passage through the polarizing neutron guide 2 and the 
filter 3 is directed along by the coil 4, then it turns round by the angle a by the coil 5. 
If the crystal does not influence the spin orientation, the polarization vector will be restored 
in the initial direction along the by the coil 8. The behavior of the neutron spin for the 
case of a = 90° is shown in Fig. [7| We use the same coordinate system (X, Y, Z) in Fig. [7] 
and Fig. |1] 

The dependence of a neutron beam intensity on the angle a after the analyzer 10 was 
studied to observe the depolarization effect for diffracted neutrons. The analyzer 10 trans- 
mits the neutrons with the polarization parallel to only. The measurements are similar 
to those, using a spin-echo technique. 

The neutron beam passed through the polycrystalline BeO filter 3 of a 120 mm thickness 
to reduce the contribution of the background reflections (see Fig-EJ to the forward diffraction 
beam. The residual contribution of them was estimated to be ~ (20 ± 10)% of the useful 
intensity of neutrons diffracted by the (110) plane. The uncertainty of this contribution 
results in a systematic error of a measured value. 
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Figure 7: The behavior of the diffracted neutron spin for the case a = 90°. 
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Figure 8: The dependence of the neutron intensity N on the angle a° between the Schwinger 
magnetic field and the vector of polarization of incident neutrons for the Bragg angle 
6» B = 84°. 
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If the neutron spins turn round by the angles ±A<^q in the crystal for the waves described 
by and ip^ 2 \ the counting rate N of the neutrons after the analyzer will be: 



N = N (l + P z ) = N (l + P (cos A4>q sin 2 a + cos 2 a)), 



(11) 



where Pz is the projection of a neutron polarization after the crystal on the direction H^. 
In the case of A0q = we will have Pz = Po and N will not depend on the angle a. The 
value of initial polarization was Po = (87 ± 3)% for neutrons with the wavelength A ~ 4.8A. 

The example of the dependence N(a) is shown in Fig. [SJ The values of polarization Pz 
are shown on the left axis of ordinates. The solid curve in Fig. [HI is a result of fitting the 
experimental points by the dependence 1)1 

As it has been noted earlier IHUHl, the angle of neutron spin rotation due to Schwinger 
interaction does not depend on the Bragg angle, and that is experimentally proved (see 
Fig. 03). 
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Figure 9: The dependence of the angle A<^g of a neutron spin rotation due to Schwinger 
interaction on the tangent of Bragg angle. A and B are two crystal positions (see Fig. 
and Fig. HJ. 

The experimentally observed result corresponds to the interplanar electric field, acting 
on a diffracted neutron, equal to 

£ ( iio) = (2.24 ± 0.05(0.20))10 8 V/cm, (12) 

The systematic error caused by uncertainty of the contribution of background reflections is 
pointed in parentheses. 
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The experimental values are in good agreement with the earlier theoretical predictions 
and confirm the opportunity to increase more than by an order of magnitude the sensitivity 
of the method to neutron EDM, using the angles of diffraction close to 90°. It is experimen- 
tally shown that the value Et determining the sensitivity of the method in our case can reach 
~ 0.2- 10 6 V s/cm, what is comparable with that of the UCN method (~ 0.6- 10 6 V s/cm)|S] 
and much more than the value obtained by Shull and Nathans (~ 0.2 ■ 10 3 V s/cm')[TT]. 

5 Experimental test of the sensitivity 

In 2002 the test experiment was carried out at ILL (Grenoble, France) pTT] . The main 
purpose of this experiment was to estimate the statistical sensitivity of the Laue diffraction 
method for the neutron EDM, using the available quartz crystal and the facility for particle 
physics with cold polarized neutrons PF1A or PF1B at the ILL high flux reactor. The 
scheme of the experiment was similar to that of the previous one (see Fig. if)). The TOF 
distance was 1 m. The experiment was carried out using the (110) plane (d=2.456A) of 
a quartz crystal with the sizes 14 x 14 x 3.5 cm 3 prepared and tested at PNPI (Gatchina, 
Russia). The mosaicity of the crystal was less than 1" over all crystal volume. 

Examples of experimental TOF spectra of the forward diffracted beam are shown in 
FigEB The positions and widths of the different reflection peaks coincide with the theoret- 
ical expectations. 

One can see from Fig^J] that the reflection peaks are located on a background. Its 
spectrum is proportional to that of the incident beam. This background is due to incoherent 
scattering of the incident neutrons inside the crystal. It can be reduced essentially by 
increasing the distance between the crystal and detector (see Fig HT) and will be negligible 
for the final geometry of the experiment. 

The measured dependencies of the time the neutrons stay inside the crystal are shown 
in Fig^| The two curves correspond to the symmetric crystal positions A and B with the 
same Bragg angles, see Fig.|2p. The difference of these two curves is due to an inaccuracy 
of the initial angular position of the crystal. The time of stay reaches 1.8 ms for the Bragg 
angle 0b ~ 88.5° which corresponds to the mean neutron velocity Vii ~ 20 m/s in the crystal. 
The incident neutron velocity was « 800 m/s (A « 5 A). 

The dependence of the (110) reflection intensity on the Bragg angle is shown in Fig 1 131 
One can see a good coincidence of the theory with the experimental data for Bragg angles 
smaller than 86°. The disagreement between theory and experiment for Bragg angles larger 
than 87° can be explained by the decrease of the area of the incident neutron beam '"seen"' 
by the crystal for large angles of diffraction. 

These measurements allow us to determine the dependence of the statistical sensitivity 
for the neutron EDM on the Bragg angle (see Fig ll4|) . This dependence has maximum for 
the Bragg angle equal to 86°. 

Using the experimentally measured values of the time of neutron passage through the 
crystal (FiglT2*)). the intensity of the diffracted beam fFig |T5fr . and the previously measured 
value of the electric field for the (110) crystallographic plane (E g = 2.2 • 10 s V/cm |3()l 132] ) 
we can estimate the statistical sensitivity of this method to measure the neutron EDM, 
using the world's highest intensity cold neutron beam PF1B and the quartz crystal. This 
sensitivity is equal to ~ 6 • 10~ 25 e • cm per day for the Bragg angle 9 b = 86°. 
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Figure 10: Spectrum of the forward diffracted neutron beam for different Bragg angles. 
Note the principal difference in the behavior of the working plane (110) and the other 
crystallographic planes (for instance, (101)). 
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Figure 11: Spectrum of the forward diffracted neutron beam for two time of flight distances. 
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Figure 12: The dependence of the time of neutron stay in the crystal on the Bragg angle. 



N, 1/s 

10 



0,1 



0,01 



80 



82 



84 



- 


■ : 


















■ \ 




Theory for 






a>=3 10 8 n/(cm 2 As) 




■ 

\ 


■ 












i 


i 




\ 

■ i 



88 90 
B 



Figure 13: Dependence of the intensity of the (110) reflection of the quartz crystal on the 
Bragg angle. The two curves correspond to the symmetrical crystal positions A and B (see 
Fig. if)). 
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Figure 14: The measured dependence of the sensitivity of the method for the neutron EDM 
on the Bragg angle. 

The comparison with the magnetic resonance method using ultracold neutrons (UCN- 
method) 011 EJ is given in Tabled 

Possible further progress of this method for the neutron EDM search experiment may 
be associated with the use of some other crystals. In principle, there are the crystals that 
can allow to increase the method sensitivity about by an order of value in comparison with 
the quartz one (see Table. 0). 
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Tabic 1: The comparison of the Laue diffraction method with the UCN one. The intensity 
for the Laue diffraction scheme is recalculated from the experimental value (see Fig. I13|) for 
the quartz crystal dimensions 3.5 x 14 x 25 cm 3 and Bragg angle 9b — 86° (it/2 — 9b — 1/15). 





UCN-method 5 


Laue diffraction 
method 


E (kV/cm) 


4.5 


2.2 -10 b 


t(s) 


130 


0.7 TO" 3 




(v =5-6 m/s) 


(«H « 50 m/s) 


Et (kV s/cm) 


585 


150 


N (neutrons/s) 


60 


1 • 10 3 


<7d e-cm per day 


6 • 10~ 2b 


6 • 10~ 2b 



Table 2: Parameters of some noncentrosymmetric crystal suitable for the EDM experiment. 
T a is a time of neutron life in the crystal (time of absorption) . 



Crystal 


Group 


hkl 


d, (A) 


E g , 


Ta, 


EgT a , 




Symmetry 






10 9 V/cm 


ms 


(kV s/cm) 


a-quartz 


32(D§) 


111 


2.236 


0.23 


1.0 


230 


(SiOa) 




110 


2.457 


0.20 




220 


Bii2Si02o 


123 


433 


1.75 


0.43 


4.0 


1720 






312 


2.72 


0.22 




880 


Bi 4 Si 3 0i2 


-43m 


242 


2.1 


0.46 


2 


920 






132 


2.75 


0.32 




640 


PbO 


Pca21 


002 


2.94 


1.04 


1 


1040 






004 


1.47 


1.0 




1000 


BcO 


6mm 


011 


2.06 


0.54 


7.0 


3700 






201 


1.13 


0.65 




4500 
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6 Neutron optics in noncentrosymmetric crystal 



Here we consider the neutron-optic effects for a neutron, moving through a noncentrosym- 
metric crystal with the energy and direction far from the Bragg ones, when the deviation 
from the exact Bragg condition reaches (10 3 — 10 5 ) Bragg widths. 
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Figure 15: Neutron movement with the different directions of the wave vectors K with the 
respect to reciprocal lattice vector g. a) |K + g| > K, neutrons are concentrated at the 
"nuclear planes" (maxima of nuclear potential), b) |K + g| < K, neutrons are concentrated 
between the "nuclear planes". These two cases correspond to neutron optics, c) The third 
case |K + g| = K corresponds to the neutron diffraction, when the both kinds of waves 
excite in the crystal. 

The essence of the phenomenon is the following. Let a neutron is moving through the 
crystal and the Bragg condition is not satisfied for any crystallographic plane. In this 
case the distribution of the neutron density IVK 1 ")! 2 in the crystal can be written using a 
perturbation theory 

|^(r)| 2 = 1 + E gfc 2 _\ cos (g r + ( 13 ) 

where E). — H 2 k 2 /2m and E/. — h 2 \h + g| 2 /2m are the energies of the neutron with the 
wave vectors k and k + g in the crystal, the values v g and 4> g are respectively the absolute 
magnitude and the phase of the (/-harmonics amplitude V g — v g expi(j} g of the neutron 
crystal interaction potential, which has a form 

V(r) = V a expi(gr) = V + ^ 2v g cos (gr + fl ), (14) 
a g 
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difference A g = Ek — Ek g describes the deviation from the Bragg condition measured in 
the energy units. One can see that the neutrons are concentrating either on the maxima or 
on the minima of the periodic potential, depending on the sign of A g (see Fig. I15ll6fl . the 
degree of this concentration being determined by the value of v g / A g . 




\K+g\<K 
\K B +g\=K B 
\K+g\>K 



Figure 16: Neutron movement with the different absolute values of the wave vectors K 
(different wave lengths). The same three cases as in the Fig. 1151 |K + g| > K, |K + g| < K 
and |K + g| = K 

The neutron interacts with the nonmagnetic crystal by the nuclear forces, therefore the 
neutrons will concentrate on the maxima (or on the minima) of the nuclear crystal potential, 
so the neutrons in the noncentrosymmetric crystal will move under electric interplanar field 
as in the diffraction case analyzed above. The sign of the electric field depends on the sign 
of the deviation parameter A ff = Ef. — Ek fl - 6 

Due to such concentration the neutron kinetic energy in the crystal Ek = h 2 k 2 /2m (in 
the second order of the perturbation theory) will be equal to 

2m 2m ^ A ~ 2m ^ A, ' 1 ' 

where fco is the wave vector of an incident neutron and k is the wave vector of the neutron 
in the crystal, K is the wave vector in the crystal with the mean nuclear potential Vo taken 
into account h 2 K 2 /2m — h 2 k 2 /2m — Vq. As follows from l|15fl . the neutron, moving through 
the crystal far from the Bragg directions, nevertheless, "feels" the crystal structure. 

6 The equality E k = E kg corresponds to exact Bragg condition for the plane system g. In this case the 
perturbation theory (valid for v g <C (E^ — i?fc s )) becomes inapplicable, and one should use the two wave 
diffraction theory. 
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For the case of non-magnetic and non-absorbing crystal the expression for V g can be 
written as |33j: 

Vg = v»e^ +iv*j<t>h^^, (16) 

where , (p^ are respectively the absolute magnitude and the phase of the g-harmonics 
amplitude of the nuclear neutron-crystal potential, ,<t> g are the absolute magnitude and 
the phase of the g-harmonics amplitude of the electric potential of crystal, fx, v are the 
magnetic moment and the velocity of neutron, c is the light speed. 

By substituting this expression into (|15fl and taking into account that for non-absorbing 
crystal V g = V2 g , we shall obtain 

tfk 2 h 2 K 2 ^«) 2 o-[E sum xv] 



where 



2m 2m ^— ' A 



^^<gsin(A0 9 ) (18) 



is a resultant electric field affecting a neutron in the crystal. Here A<j) g = (f) 1 ^ — (j>^ is the 
phase shift between g-harmonics of nuclear and electric potentials of the crystal. 

We should note that the neutron refraction index n in this case will depend on the 
neutron spin direction 

n 2 = k 2 /kl = nl - An 2 , - An*, (19) 
where Uq = K 2 /k 2 is a square of a mean refraction index for neutron in the crystal, 

An 2 d = (2m/h 2 k 2 )J2(v g V ) 2 /A g (20) 

9 

is a small diffraction correction to the square of mean refraction index, 

An 2 s = (2m/h 2 k 2 )([MT[E sum x v])/c (21) 

is a spin dependent correction arisen due to interference of nuclear and electric amplitudes. 

For the centrosymmetric crystal A<f> g = and therefore E SUI „ = 0. In the noncentrosym- 
metric crystals there are the crystallographic planes, for which A(j> g ^ and so the electric 
field acting on a neutron will be nonzero E snm =/= 0. 

Therefore a spin dependence arises for the neutron-crystal interaction, which leads to 
different values of a neutron wave vector in the crystal for two opposite spin orientations 
that in turn leads to neutron spin rotation around the direction of Schwinger magnetic field 

H Ln = [Esum X V]/C. 

The rotation angle for the crystal length L will be equal 

2^[E. um xv]L 

h c v 

One should add an imaginary part into the nuclear crystal potential to describe the 
absorbing crystal: 

Vg = v N e i*» + iv N' e i<t>? + iv E e i4>f^^A_ (23) 
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Here v 



N' ±N r 



are the amplitude and phase of g-harmonics of the imaginary part of the 



nuclear potential. 

The value of kinetic energy in the crystal becomes equal to: 



where 



2 1,2 



2m 



h 2 K 2 
2m 



v {g) -i(v{ + v( g) )-» 



c[(E sum + lEi' sum ) X v] 



Vi 



{a) 



E 



v «) 2 -«') 2 

g y 
2v?vf cos« - <f>f) 



(24) 

(25) 
(26) 



E' 



E 



2v N ' 



(27) 



The estimations give that the values of the diffraction corrections to a mean potential 
for a-quartz crystal are + iV/ g \ ~ 10~ 3 (Vb + iVq), 

H<r[(Ei sum + i~E' sum ) x v]/c ~ 10 _6 (Vo + iVq) for the wide range of the incident neutron 
wavelengths and sharply increase near the Bragg conditions. We should note also that in 
spite of a smallness the last correction leads to relatively large and observable effects due to 
its spin dependence. 

7 Observation of the neutron spin rotation 

The experiment was carried out at the PNPI WWR-M reactor. Scheme of the experiment 
is shown in Fig.[T3 see [H EH EU . 




Figure 17: Scheme of the experiment. 1 is a polarizer; 2 is a 7r/2 coil to turn spin around X 
axis; 3 is the a-quartz single crystal with the sizes 14 x 14 x 3.5 cm 3 ; 4 is a ±7r/2 coil to turn 
around Y axis; 5 is an analyzer. is the guiding magnetic field; O is an axis (parallel to 
Z direction) of a crystal rotation; and P f are the polarizations of neutron beam before 
the crystal and after it. 
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Initially the neutron spin was directed along the neutron velocity (axis Y). The X- 
component of the polarization was measured after neutron passage through the crystal. 
This component should be equal to zero, if the spin rotation effect is absent. Time of flight 
technique was used for measuring the spectral dependence of polarization. The crystal 
was overturned around Z axis to eliminate the false effect due to nonzero value of the X- 
component of polarization for real setup. The effect changes its sign due to a change of the 
sign of the electric field in this case. 

The measurement was carried out using the a-quartz crystal with the dimensions 14 x 
14 x 3.5 cm 3 . The crystal orientation was determined by the angle (3 between the neutron 
velocity (Y axis) and the [110] crystal axis. 

The theoretical dependence of the angle of neutron spin rotation A(p s on the value and 
direction of neutron wavevector is shown in Fig. 




Figure 18: Theoretical dependence of A(^ s on the wavevector and direction of neutron for 
the a-quartz crystal. 

The experiment was carried out for two crystal positions with (3 = 90° and 30°. The 
results are shown in Figs. [5D| Two plots at the figures correspond to different energy 
resolution of the experiment (AA/A = 5 • 10~ 2 and = 2 • 10~ 2 ). The solid curves reproduce 
the theoretical dependence 122|) averaged over energy resolution. The dotted lines indicate 
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the positions of the crystallographic planes with nonzero value of A(j> g (see ((THJ) . One can 
see a good agreement between theoretical and experimental results. The resultant electric 
field |E sum | is shown on the right ordinates axis. Its value is |E sum | « (1 - 10) • 10 4 V/cm 
for any point of spectrum, and so Aip s may reach ±2 • 10 -4 rad/cm. 



A<p , rad/cm 



E V/cm 

sum ' 




Figure 19: Energy dependence of Aip s for (3 = 90°. Solid curves are the theoretical depen- 
dence (|22|l after averaging over the energy resolution of the experiment. 



8 Spin rotation effect for Bragg reflected neutrons from 
the deformed crystal part 

In the previous chapter we have considered the effect of neutron spin rotation for large 
deviations from exact Bragg condition (~ 10 3 — 10 4 Bragg width). The measured effect has 
coincided with the theoretical one but has a small value |.'S5| . Here we consider the situation 
when we can tune the sign and value of the deviation from the exact Bragg condition within 
a few Bragg width. In this case the neutron spin rotation angle essentially increases due to 
increasing the electric field affected the neutron. 

Let's consider the symmetric Bragg diffraction scheme with the Bragg angles close to 
the right one. Neutrons fall on the crystal in the given direction with the energy close to 
the Bragg one for the crystallographic plane g. Deviation from the exact Bragg condition is 
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A<p [rad/cml E [V/cml 

s sum J 




Figure 20: Energy dependence of Aip s for (3 = 30°. 

described by the parameter A g — Ek — E^ , where Ek — h 2 k 2 /2m and Ek. = fi. 2 |k + g| 2 /2m 
are the energies of a neutron in the states \k) and \k + g) respectively. 

In this case the neutron wave function inside the crystal in the first order of perturbation 
theory can be written |33| 

ip(r) =e-' lkr + a-e- l(k+s '> r , (28) 

where 

3 m _ \v a \_ 

Ek — Ek g A g 

Here V g is g-harmonic of interaction potential of neutron with crystal, 
consider the case a <C 1, so we can use the perturbation theory. 
The electric field affected the diffracted neutron will be equal to 

E = 2E g • a, (30) 

where E g is the interplanar electric field for the exact Bragg condition. 

One can see that the sign and value of the electric field (|30|l are determined by the 
sign and value of deviation A g from the exact Bragg condition, therefore to have the given 
electric field and so the effect of neutron spin rotation we should select from the whole beam 
the neutrons with the corresponding deviation parameter A g . 



(29) 

For simplicity we 



2G 



The presence of the electric field will lead to an appearance of the Schwinger magnetic 
field 

U s = l/c[E x V||]. (31) 
The neutron spin will rotate around the H5 by the angle 

Vs = —z , (32) 

nv± 

L c is the crystal thickness, V|| and vj_ are the components of neutron velocity parallel and 
perpendicular to the crystallographic plane correspondingly. 

The main idea of the experiment is the following. We use a small controlled variation 
of the interplanar distance Ad (caused by heating, for example) near the exit crystal edge. 
So some part of neutrons passed through the crystal will reflect from this small crystal 
part. These back diffracted neutrons have the deviation parameter for the main part of 
crystal determined by Ad and so they propagate under corresponding electric field in both 
directions there and back. Thermal deformation of the crystal edge is used to create such 
variation of the interplanar distance. 



v, - -i T=T ±AT 
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Figure 21: Two crystals are in parallel position. Neutron reflected by the small crystal 
pass twice through the large crystal. The deviation parameter A g for the large crystal is 
determined by the temperature difference AT. 

We can use also two separate crystals in parallel position for this purpose (see Fig l21fl . 
One can heat (or cool) the second small crystal. The part of neutrons passed through the 
first crystal with the corresponding Bragg wave length will reflect by the second crystal with 
the given deviation parameter for the first (large) crystal. This deviation parameter will 
directly depend on the temperature difference between crystals. 

Value of the wave length Bragg width for (110) quartz plane (d = 2.45A) is AAs/A ~ 
10~ 5 . To shift the wave length of the reflected neutrons by the one Bragg width we should 
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heat (or cool) the second crystal to have the same value of Ad/d. Linear coefficient of the 
thermal expansion for quartz crystal is AL/L as 10~ 5 per degree. Therefore, the deviation 
±AAb corresponds to difference of the crystal temperatures AT w ±1°. We note that the 
different signs of this temperature difference will correspond to different signs of the electric 
field acting on the neutron. 




(110) Peltier 
element 

Figure 22: Passage of the neutron through the crystal. Presence of the interplanar distance 
gradient result in forming the reflex near the back face of crystal. 

Scheme of the neutron behavior in the crystal is shown in Fig l22l Two samples of quartz 
crystal has been used in this experiment with the thicknesses along X axis L c — 14 and 27 
cm. The Peltier element has been attached to the back face of the crystal. That allows to 
create the temperature gradient in the crystal along the neutron trajectory. So the Bragg 
condition will vary along the neutron trajectory and different parts of crystal reflect the 
neutrons with the different A. Therefore, the reflected beam will contain not only the reflex 
from the entrance crystal face (corresponding to Bragg condition for d) but also the reflection 
from the back exit face (corresponding d ± Ad) that twice pass through the crystal there 
and back. Moreover, the value of the deviation parameter A g for this reflection is directly 
depend on the value of temperature gradient. In the case of higher temperature of the back 
crystal face the neutron with Ek — Ek g > will be reflected, while in the case of its lower 
temperature the neutron with Ek — Ek < will be reflected. 
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Examples of the time of flight spectra of the reflected neutrons for the Bragg angle ~ 90° 
are shown in Fig[33] One can see a formation of the reflex from the back crystal surface 
and increasing its intensity with the rise of the temperature gradient. 




0,000 



1.0 1.5 
TOF, ms 



Figure 23: Dependence of the time of flight (TOF) spectra of the neutron reflected by 
the (110) plane of quartz on the temperature gradient applied to the crystal. Bragg angle 
~ 90°, L c = 27 cm. One can see the reflexes from the front surface and from the back part 
of crystal. 

The scheme of the experiment on the observation of neutron spin rotation is similar to 
that described in |32| . 

To observe the effect of neutron spin rotation due to Scwinger interaction it is necessary 
to turn the crystal in a position, for which Bragg angle is different from 90°, because in the 
case of Bragg diffraction the Schwinger effect disappears for 90° Bragg angle: 

4E /^cV|| 4E^L C 

tp s = = T—ctg(6 B ) — ► (33) 

cnv± cn e B ^ir/2 

The experiment on the observation of neutron spin rotation was carried out with the 
L c = 14 cm crystal thickness and Bragg angle ~ 86°. 



29 



The dependence of the angle of neutron spin rotation around H5 on the value of tem- 
perature gradient is shown in Figl24l 
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Figure 24: The dependence of the angle of neutron spin rotation due to Schwinger interaction 
on the value of temperature gradient. Two upper figure corresponds to two crystal positions 
differing by the angle 180°. One can see a good coincidence of the theoretical dependence 
(solid curve in the bottom plot) with the experimental points. 

We can change the sign of the effect by turn the crystal by the 180° around H5. One 
can see that the experiment confirms that such a crystal rotation indeed change the sigh of 
the observed effect. On the right axis the effective electric field that is necessary to get the 
corresponding spin rotation effect is shown. One can see that the value of the electric field 
reaches ~ 1.3 • 10 s V/cm, that is only 1.5 times less than in the Laue diffraction case for 
exact Bragg condition EOj • 
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9 Conclusion 



The first experimental study of some new phenomena for the neutron diffraction and optics 
in the noncentrosymmetric crystal was carried out, using the pilot set-up created for a search 
for the neutron EDM by the crystal-diffraction method. 

For the first time the neutron dynamical Laue diffraction for the Bragg angles close to 
the right one (up to 87°) was studied, using the forward diffraction beam and the thick 
(~3.5 cm) crystal. 

The effect of the essential time delay of diffracting neutrons inside the crystal for Bragg 
angles close to 90° was experimentally observed. For (llO)-plane of a-quartz and 9b = 88.5° 
we have obtained tl ~ 1.8 ms that corresponds to uii « 20 m/s, while v = 808 m/s. 

The predicted earlier phenomenon of the neutron beam depolarization was first exper- 
imentally observed for the case of Laue diffraction in the noncentrosymmetric a-quartz 
crystal. It is experimentally proved that the interplanar electric field, affecting the neutron 
in the crystal, maintains its value up to Bragg angles equal to 87°. 

It is shown experimentally that the value Et determining the sensitivity of the method 
in our case can reach ~ 2 • 10 5 V s/cm that is comparable with that for the UCN method 
(~ 6 • 10 5 V s/cm) |i>| and much more than the value obtained by Shull and Nathans (~ 
2 • 10 2 V s/cm) HU and also than that of 36 (~ 1.2 • 10 3 V s/cm) 7 . 

These results give the opportunity to try the Laue diffraction method for a neutron 
EDM search. The statistical sensitivity of the method was estimated experimentally to be 
cr(D) «6x 10~ 25 e • cm per day for the PF1B beam of the ILL reactor and available quartz 
crystal. The measured intensities and the neutron time of stay in the quartz crystal coincide 
well with the theoretical predictions. The use of the other crystals may allow to improve the 
sensitivity of the method by about one order of magnitude and to reach ~ a few 10 _26 e- cm 
per day. 

The effect of a spin rotation due to Schwinger interaction of the magnetic moment of 
moving neutron with an interplanar electric field of the noncentrosymmetric crystal was 
experimentally observed for neutron, passing through the crystal far from the Bragg condi- 
tions. 

The energy dependence of a spin rotation angle was measured for two crystal orientations. 
The dependence has a pronounced resonance form. The direction of a spin rotation changes 
its sign when the energy passes through the Bragg resonance value. So a new kind of 
neutronography arises for noncentrosymmetric crystals, which allows to "see" and study 
the crystallographic planes with nonzero electric interplanar fields. 

For a-quartz crystal the value of the spin rotation angle can reach ±(1 — 2) • 10 -4 rad/cm 
that corresponds to the value of resultant electric field equal to ±(0.5 — 1) • 10 5 V/cm. 

We note also that the presence of the terms <|26|) and (|27fl in equation l|24|) should result 
in the dependence of a neutron absorption on the direction and the value of a neutron wave 
vector as well as on the neutron spin orientation. 

The effect of spin rotation in a noncentrosymmetric quartz crystal for neutrons Bragg 
reflected by the deformed part of crystal was first observed. This effect is caused by the 

7 One can consider the measurements of the depolarization as a first preliminary and rough measurements 
of a neutron EDM, which gives the result D < 10~ 22 e-cm that is some better than the old Shull's and 
Nathans's result 1111 . The result obtained in magnetic resonance method using cold neutrons 1361 was 
D < 3 ■ 10" 24 e-cm. 
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Schwinger interaction and depends on a deformation degree of crystal near its back surface. 
For the quartz crystal the effective electric field affected the neutron during the time of its 
staying inside the crystal can reach ~ 1.3 • 10 8 V/cm. Simple estimation has shown that 
in our case the depth of neutron penetration into the crystal and so the time of neutron 
interaction with the electric field can be about four or even five orders more than in the well 
known Shull and Nathans experiment for the neutron EDM search |11|. 

In addition, the requirements to the crystal perfection are relatively low for this scheme. 
For the case 75 << w m the effective electric field affected the neutron depends on an 
effective crystal mosaicity w m as E — Eq(jb /w m ), where 75 is the angular Bragg width, 
but the reflex intensity increase as / = Io{wm/jB), therefore the sensitivity to measure the 
neutron EDM will be reduced only by a factor yJw m /"fB, that give us a hope that such 
a scheme can be applied to search the T-odd part of neutron-nuclei interaction [2J using 
neutrons with energies near the P wave resonance one. 
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